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ABSTRACT 

The evidence that variations in the sun itself are 
responsible for climate changes is substantial and hence it 
is frequently suggested that the sun's near-ultraviolet 
radiation is significantly variable. Whether or not that 
inference is true, the evidence is strong that the sun's 
magnetically associated phenomena have long-term variations, 
aside from the familiar 11-year cycle. A long period of low 
solar activity and enhanced cosmic-ray production of in 
the stratosphere will reduce the ozone abundance. 

A perturbation analysis, allowing for temperature and 
opacity feedbacks, is developed to calculate depletions in 
the abundance and reductions of stratospheric solar 
heating that result from increases in NO^ concentration. A 
pair of "perturbation coefficients" give the reduction in 
0^ and temperature through the stratosphere for a specified 
NO increase. This type of analysis illustrates the tendency 
for various levels to self-heal when a perturbation occurs. 

Physical arguments appear to indicate that the expected 
sign of the climatic effect is correct, with colder surface 
temperatures produced by reduced magnetic shielding. In 
addition, four qualitative reasons are suggested for thinking 
significant ozone reductions by cosmic-ray influxes will 
lead to an increased terrestrial albedo from stratospheric 
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ccndonsation. In this view long-term ('vlO years) climatic 

changes have resulted from secular geomagnetic variations 

2 

while shorter ('\,10 years) excursions are related to changes 
in solar activity. 
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1. Introduction 

Climatic changes have been attributed to internal 
exchanges of energy (among atmosphere, oceans, ice caps) as 
well as to external modulations, and even the relative 
imx^ortance of the two classes of mechanism is not settled 
(Schneider and Dickinson, 1974). Nevertheless, a scholarly 
historical study by Eddy (1976) provides convincing evidence 
that the coldest part of the "little ice age", between about 
1650 and 1715, coincided with a near absence of solar activity 
(the "llaunder minimum"), 

Schneider and Mass (1975) have obtained a striking 
agreement between the global surface temperatures since 1600 
and temperatures computed with a simple formula based on 
radiative equilibrium. This approach relates the equilibrium 
temx>erature to the solar constant and then connects the solar 
constant to the sunspot number by the Kondratyev-Nikolsky 
(1970) empirical formula, which has no theoretical justification 
but was thought to describe 11-year variations in the solar 
constant. This sunspot formula implies an increase in solar 
flux up to 2.5 percent as the Wolf sunspot number increases 
from 0 to 80 (indicating moderate activity) and then slowly 
decreases again to 1.9 percent enhancement for VJolf numbers 
of 200 (high activi*y) . As Schneider and Mass emi^hasize, it 
is not important whether the formula gives a correct relation 
between the solar constant and sunspots; it may merely serve 
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as a link in a convenient relationship between surface 

teiiiperature end solar activity. The fact that the link, goes 

through solar flux as the physical connection does not mean 

it is so. Any effect, such as solar particle emission, solar 

magnetic fields, or atmospheric transparency (or, as we shall 

argue below, all three) may just as well be the physical 

connection between sunspot and climate. 

Suess (1965, 1968, 1971, 1974) has shov?n, by two 

independent correlations, that climate is associated 

v/ith solar activity over the past 6000-3000 years. First, 

14 

production of radioactive C (by galactic cosmic rays) 

fluctuates due to changes in solar magnetic activity. Second, 
14 

the C level seems to fluctuate v;ith climate, in the sense 

that prolonged high solar activity leads to warm winters. 

Long-term changes (of the order of 10^ years) in the 
14 

C level are likely due to a change in the geomagnetic dipole 
moment (Ramaty, 1965; Suess, 1965). For still longer periods 
it seems reasonable to conclude that faunal extinctions that 
are correlated with reversals of the earth's magnetic field 
over the past 2.5 million years (Hays, 1971; Opdyke , 1972) 
are due to climatic cltanges induced by the increased cosmic-ray 
deposition in the middle atmosphere (Harrison, 1968), rather 
than by direct cosmic-ray effects at the earth's surface 
(Uffen, • 


Nevertheless, it has not been clear what physical 
mechanism would relate cosmic rays to the Earth's climate. 

In the meantime a rather different line of study has indicated 
a modulating influence of the sun on the stratosi^here, not 
through variable solar radiation but through the sun's 11-year 
cycle of solar-wind activity. 

Warneck (1972) and Brasseur and Nicolet (1973) noted 
that cosmic-ray ionization in the stratosphere produced free 
N and N^, which in turn react with oxygen to form HO^. 

Although this source of NOj^ (which catalytically destroys 
ozone) is relatively minor, Ruderman and Chamberlain (1975) 
showed that the observationally known 11-year variation in 
cosmic-ray ionization over the polar cap implied a predictable 
11-year variation of ozone, both over the polar cap and, with 
some months' time lag, at lower latitudes as well. Analysis 
of the available data (Angell and Korshover, 1973, 1975) 
strongly supports the existence of an 11-year ozone cycle with 
approximately the predicted phase lag from sunspot maximum. 

The observed amplitude seems to be stronger than was expected 
from simple ionization and subsequent production of NO . 

Hence there is reason to think additional chemical reactions 
occur - perhaps involving negative ions - that could be 
expecially important in destroying ozone over the polar caps 
(Ruderman, Foley, and Chanil^erlain, 1976). 
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The coitibination of the liistorical relationship between 
solar activity and climate with a chemical link between 
low-energy cosmic rays and the ozone balance may t>e the 
crucial clue; It suggests that heliomagnetic and geomagnetic 
control of cosmic**ray fluxes and ozone destruction in the 
stratosphere is a principal factor influencing climate. 

Alternative or additional causal mechanisms are 
nevertheless possible. Ruderman (1974) has suggested that 
faunal extinctions could be induced by explosions of nearby 
superjiovae. The dramatically increased x-ray or cosmic 
radiation at tlie earth could reduce the protective ozone 
shield, subjecting life forms (especially DIIA) to increased 
ultraviolet exposure. Life v;ould be especially vulnerable to 
radiation from giant solar flares during this period. 

Reid et al . (1976) have proposed that faunal extinctions 
occur during times of geomagnetic field reversals because 
the stratosphere is then more susceptible to the occa“i ‘'al 
solar cosmic rays. TiS with Ruderman' s supernova, the particle 
bomljiardment creates stratospheric WOjj, which depletes ozone, 
raaj'.ing surface life susceptible to solar ultraviolet radiation. 

Whether or not life has been directly subjected to 
hazardous increases in ultraviolet radiation, long-period 
changes in tlie heliomagnetic or geomagnetic field seem certain 
to affect the ozone abundance in the stratosphere. In the 
next section we develop a perturbation treatment to estimate 

REPROBIjCr OF tdh, 

UKAO iiih' ^ t - 
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the redistribution with hei nd the magnitude of total 
ozone decrease accompanying a rise in NO^.. The concluding 
section examines the likely influence on climate, including 
the difficulty posed by certain "feedback" mechanisms. We 
propose that the stratosphere itself is the primary cause 
climate changes through increased water condensation. 

2, Perturbation Treatment of Stratospheric 0^ Abundance 

Consider an oxygen-nitrogen atmosphere. The principal 
reactions with rate coefficients, k, and dissociation rates 
at 80 nm, J(80), are: 


HO + O 3 

* HO 2 + 0 ^, ki = 2 X 10 ”^^ g-1225/T 

cm /s , 

( 1 ) 

NO 2 + 0 

-12 3 

NO + O 2 , k 2 = 9.12 X 10 cm /s , 


( 2 ) 

HO 2 + h^ 

“ NO +0, J3(80) = 4.1 X 10~^ s"^. 


(3) 

O 3 + 0 

^ -2140/T 

O 2 + O 2 , k 4 «= 1 . 3 X 10 e 

3 , 

cm /s , 

(4) 

0 + O 2 

+ M O 3 + M, kg = 1.05 X 10 g520/T ^ 

(5) 

O 2 + h ’ 

-9 -1 

-* 0 + 0 , Jg (80) = 4 X 10 s , 


( 6 ) 

O 3 + h.^ 

-3 -1 

-* 0 + O 2 , J 7 (80) = 4 . 8 X 10 s 


.7) 


The rates of reaction and photolysis are those quoted by 
McElroy et al . (1974). (Photodissociation (3) of UO 2 is 
negligible as a source of 0 atoms.) 

In the Earth's stratosphere, reactions involving hydrogen 
(OH, HO 2 / etc.) are important at the liigher altitudes. Also, 


is an inportant component of the total IJOj^ in the lower 
stratosphere, but adjustments in relative abundances of NO 
and NO 2 to the total UOjj can be made as an afterthought. 

Given the [0] , IO 3 I concentrations the [NO] ratio in 
lihotochemical equilibrium is 


r , . t_N 0 ] 

- [NO] [N02] 


[0] k2+ J3 

[^]^ kY+ M ^2+^3 


( 8 ) 


and, similarly, r 2 = 1 - r^. If r^ is known, the [0] and 
[ 02 ] concentrations are solutions of quadratic equations. 
Define 

A = k. (k, [NO] + 2J^) , (9) 

B = k2[N02l (kj^tNO] + J^) + kj^lNO] k^ [02llM] 

+ 2k^ [02] , (10) 


and 


C ^ - 2 k^ [O2] tM] [O2] 


( 11 ) 


Tlien v;e have 


IO 3 ) = 


(B - 4AC) - -B 
2 A 


( 12 ) 


J 


Similarly, with 


A' e k^(k2tN02] + 2 k3t02][M]) / 


U3) 


B' = k2 tN02] (kj^tNOl + J^) 

+ k^[HO] k3t02)lM] - 2 k^ Jet02l , 


and 


(14) 


C' = -(k^tNOj + J^) 2IO2) Jg , 


(15) 


we obtain 


io] 


(B 


.2 


4 A‘C') ' -B' 
~Th^ 


(16) 


Equations (12) and (16) may be solved by iteration with 
equation (8) to obtain equilibrium abundances of [NO], [N 02 ), 
[OJ , ana [O^] , when the total [HO + NO 2 ] is specified. We 
shall call this the "equilibrium solution". 

Wo*,.' let the total be varied by a relative amount 


i [NO 4 NO 2 ] 
[NO 4 N02r 


(17) 


V?e doflne the intermediate quantities 


F 


2C _ 
G^-GB 



/ 


( 18 ) 


where 


G = (B^ - 4 AC) " . 


(19) 


Then if the relative change in [NOj and [NO2 ] were known we 
would have 


6[OJ 


( 20 ) 


and 

- F' aA" - aB'/G' - aC'2A7(G'^-G'B') , (21) 

where F' , G' are defined analogously to F, G. 

In the sane way v\’e may write the perturbations 6 [WO] , 
a[n 02 J in terms of ^[0], ^[0^3 from equation (8), where 


1[M0] _ . . ‘^^l 

IGOJ r^ * 


^2 a to) ^2 ^^^3^ 
j 10] ■ (O 3 ] 


( 22 ) 


f 
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and where 


j - 1 + J2/k2l‘r'j . 


Similarly we hav^e 


(23) 


d[0] 

ii^r = ^ “ T lor ^ leg 


(24) 


Equatioru^ (20) , (21) , (23) . and (24) form a complete set of 
equations tliat fix the variations of lOj , {O 3 ] , [NO] , and 
IHO^J once 6 of equation (17) is specified. 

Ue now define 

ao^] tM] k. + tNO^l k^) )t. thOJ 

f = — £ i — _ 4 l^.gNOl F, (25) 

(J-, + [NO] )^,) )c_[NO,] 

q = , (26) 


([O^jtM] I<5 + [NO^] k^) lc^[NO] 


4 [O^ I k^ [NO] A' 


G 


G'B 


(27) 


( 


u 


and 


(J^ + tNOj kj^) k2tW02] 

+ k2 INO^] F' . ( 28 ) 

The solutions are 


*^1 

[O3I 


q£ ^ ~ fq ^ + j ( 

_r2f"” q'+ j(l+r^g - r^f) ^ 



( 29 ) 


“[OT 


j(gf'~ £q ") + j (£ ^■^ g*) 
r2f‘- t^g'+ j (l+rj^g-r2f) _ 


( 30 ) 


and 


1^1 



r2(jf + jg - f‘- g') 

r2f^“ rj^g'+ jd + r^g-r^f) 


p ^ A > (31) 


where the identities define the p-jirturbation coefficients 
)['l, . 3 , and P^. Then a [HO] and a [HO 2] are given by (22) and 

( 24 ) . 

There are nov. tv.»o impiM*tant feedljack mechanisrus to 
consider. ( 1 ) A change in [O^] and [1J02] alters the solar 
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heating, which altei*s the temx>ei'aturo, which changoa the rate 
coefficients. (2) A change in [O^J at iieight x alters the 
opacity (principally for 0^ jjliotodissociation) , which changes 
the photolysis rate J^. 

Of course, if tlie ternperature chajures, the rate of 
cooling also changes, but for small xi»erturbations the modified 
equilibrium temperature may be estiwated from Dickinson's 
(19 7 3) perturbation analysis. Let 'i,, (^K/day) be the heating 
rate per molecule due to 0^ absorivtion of sunlight, the 
same for N02^ and a (day Dickinson's cooling rate for small 
temperature changes. The rate of teiiipcrature response to 
departure- from equilibrium values is 

riT 

= £.(. 1 ,( 03 ]) + - a AT. (32) 


The heating rate ii. suffers from the overhead opacity and can 
be tabulated as a function of the integrated overhead ocone, 
1 ^, for a given model. When the pertuid^ed system again 
reaches a steady state the temperature change is 


.. 1N0,]A r, 

1 '^' " ““^“a “ ■ rj "1^ 


iO^] 

~ 


t'v-- * -'"o' 


13 31 


where 




d"t 

K* 

u 0 


(34) 


is an effective attenuation cross section for the heating 
sunlight. S:»nce will be the opposite sign of 6, the 
term in brackets represents a tendency for self-healing; A 
l^ositive 6 leads to less heating in the upper levels, but 
more sunlight penetrates to lower levels, where the diminished 
heating starts to recover. The [W02l term also tends to 
compensate the [O^] absox'ption, especially at lower levels, 
because the radiation it absorbs is essentially unattenuated 
(Ramanathan, et a1 . , 1976 ; Luther, 1976). 

Horeover, the change from (33) iiss a negative 

feedback on (r^cElroy, et al . 1974). A positive 6 gives 

an ii itial + 1 i''i equation (29) and -4 [0^ ] and -IT. Then k^ 
and k^ (destroying 0^) are decreased and k^ (producing 0^) is 
increased. The second iteration of a perturbation thus tends 
to recover from the initial change. Blake and Lindzen (1973) 
point out that this effect also accelerates thermal relaxation 
compared with the rate, a, derived from infrared cooling alone. 

Attenuation of sunlight affects directly; an effective 
attenuation cross section, analogous to equation (34), can be 
similarly defined and is numerically similar since the same 
radiation is involved. Once the first order T and *^7 
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found for a given perturbation# a revised equilibrium 
abundance, {^ 2^2 computed with equations (8), (12), and ( 16 ) 
with these first-order corrections. Comparison with the first 
correction, [0^]^, gives the response of ozone to the effects 
of wT and alone. One can then readily estimate the 
convergent values, as follows: 

Let the converged solution of LtO^l be expressed as 


^^^3^ _ 4T _ , ^ ,iT 

[O3] ” ’ [ 0 ^] “ T "*'1'’' T 


( 35 ) 


where is the "initial" perturbation given by ( 29 ) , B is 
a constant, and LT is the converged, final temperature 
change. From equation ( 33 ) the temperature change, 
is nearly proi^ortional to be written generally 

as 

«T ^103] 

T " '■ • (36) 

The first-order iteration is tlien 



11103: 


[03I 


/ 


T 


( 37 ) 


which yields i, since and known. The 

equilibrium abundance, 10212 / with the revised T and gives 



103] 



(38) 


which yields Putting equation (36) into (35) we have the 
converged solution, 


i-i (03] 

'lojr " r+TT 




(39) 


The temperature change may now be conveniently related 
to 'j by a perturbation coefficient, x , obtained from equations 
(36) and (39) : 

^ = -(-rfi-)6H - t i. (40) 

Finally an equilibrium calculation may lie performed to test 
the convergence of the various perturbations that follow from 
a change, 6 , in [NO ]. 

Another feedback effect noted by HcElr .,y et al . (19 74) 

is the level of production of 0 (^^); if O 3 is reduced, the 
natural production of from W 2 O + 0 (^£) is increased at low 

altitudes and decreased at higher ones. The overall effect 
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is small and can be ignored anyway if we regard i as the net 
change in [NO„] . 

X 

The input data and the results of this analysis are 
explained in Figures 1 and 2. A positive and t mean that 
both [Oo] and T decrease as [NO ] increases. Unfortvmately 
with the present approach the temperature decrease cannot be 
followed all the way to the tropopause, because at large 
optical thicknesses the concept of a local cooling rate 
balancing the local heating by sunlight no longer governs 
t’»e temperature. Radiative exchange among different parts 
of the atmosphere and the ground becomes dominant. The 
I^rincipal effect of decreasing 0^ in the. lower stratosphere 
is to decrease the temperature because less radiation from 
the ground is absorbed by 0^ in the 9.6um region (Ramanathan, 
et al . , 1976) . 

3. Discussion: The Effect of Stratospheric Ozone on Climate 

The usual (but not unanimous) conclusion from studies 
with radiative-convective models (e.g., Manabe and 
Wetherald, 1967; Ramanathan, et al . , 1976) is that lowered 
stratospheric ozone will lower surface temperatues . 
Nevertheless, the feedback processes are numerous and 
complex. One of the more important aspects of ozone is its 
dynamical transport and relatively inert chemistry in the 
lower stratosphere. Dickinson (1974) has cautioned that we 
know very little about how climatic changes would feed back 


on stratospheric dynamics, which must be coupled to 
circulation in the tropospihere. Evidently an increased 
surface heating by low lying absorbing particles could even 
reverse (Reck, 1976) the net cooling calculated for no 
absorbing particles and an ozone depletion (e.g. , 

Ramanathan, et al . , 1976), 

Thus one way in which reductions in ozone could 
initiate colder climatic periods is through radiative- 
convective impact of the stratosphere on the surface. It 
appears that the increased solar flux reaching the surface, 
with an ozone reduction, is more than compensated by surface 
cooling through reduction of the greenhouse effect at the 
9.6 um vibration - rotation band (Ramanathan et al . , 1976). 

A second, and primarily stratospheric, possibility 
which seems to deserve rather more attention is that a 
decrease in ozone will lead to haze or cloud formation in 
the lower stratosphere, significantly raising the Earth's 
albedo. There are several items that seem to bear on this 
possiblity . 

(1) The H^O mixing level in the lower stratosphere 
— 6 

(around 5 x 10 by volume) is close to the saturation vapor 
pressure at the tropical tropopause (c. 205 °k at 16 km). On 
the pressure scale the 0 is a factor of about 5 too low 
for condensation to occur; but on the temperature scale the 
tropic tropopause averages only about 10 above the 
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saturation point. This coincidence has led to the conunon 
belief {Harteck and Jensen, 1948; Mastenbrook, 1971) that 
the ’’cold trap" of the tropopause sets the upper limit to 
the mixing ratio there and above. Even under present 
conditions stratospheric clouds are not unknown. Newell 
(1972) believes that the stratospheric mother-of-pearl 
clouds occasionally seen over Iceland in the winter arise 
from moist air systematically forced upward in its zonal 
flow. 

(2) Not only does a decrease in ozone decrease the 
stratospheric temperature, but the temperature profile near 
the tropopause is altered. Manabe and S trickier (1964) and 
Reck (1976) find the tropopause to be erased by a complete 
(presumably catastrophic) removal of 0^ , independent of 
v;ater clouds or aerosols in the models. A change in the 
dynamical interaction of the troposphere and stratosphere 
could be most important to the moisture content and 
condensation in the stratosphere. 

(3) Condensation nuclei are, of course, important. It 
would be useful to understand how aerosol chemistry would 
change in the stratosphere with an ozone reduction. In 
addition, if the 0^ reduction is related to increased 
cosmic-ray ionization reaching the lower stratosphere, as 
proposed in this paper, the temperature decrease in 
conjunction with an increased ionization rate could be 
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responsibliSi for increasing the albedo. The ion densities in 

4 -3 

the stratosphere are the order of 10 cm (Ruderman, et al . 
1976) , with both positxve and negative ions being hydrated 
with several water molecules. Presumably larger initial 
water clusters are needed for continued growth. The role of 
these ions in attaching to aerosols before they are 
electrically neutralized (with lifetimes of about 300 sec) 
is not well understood. 

(4) Long-term feedbacks involving the oceans, ice-caps, 
and circulation of the troposphere are not easy to predict, 
and natural oscillations within the system are conceivable. 
But there is a first-order feedback that acts in the 
positive or reinforcing direction and is likely to be 
exceptionally important. A small decrease in stratospheric 
ozone probably decreases the ground temperature, whether by 
condensation of stratospheric moisture or merely by opening 
the 9.6um cooling window. The cooler ground radiates less 
energy to the region just above the tropopause , which is 
heated through infrared absorption by ozone (as well as CO 2 
and H 2 O) and not so much by ultraviolet absorption of 
sunlight. This is the crucial altitude for increased 
condensation, and the diminished tropopause teiuperature will 
lead to further condensation. And so on. 

One problem, however, that does not seem to be settled 
is whether condensation of stratospheric water would cause 
enough heating by Infrared absorption to offset the 
increased albedo to sunlight. The formation of high cirrus 




iSi 


20 


clouds in the troposphere probably does heat not only the 
tropopause but the ground (Manabe and Strickler, 1964). But 
it is not clear to what extent condensation in a nearly 
saturated stratosphere would modify the warming effects 
there compared with the warming caused by the low (absolute) 
humidity already present. 

I am grateful to Ste’^en Wofsy and Michael McElroy for 
supplying me with the detailed data that were used in 
developing their 1974 model of the stratosphere, which 
served as the starting point for the perturbation analysis 
of section 2. The research reported in this paper was 
supported by the Atmospheric Sciences Section of the 
National Science Foundation and the National Aeronautics and 
Space Administration. 
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Legends 

Fig. 1. Adopted height distribution of odd nitrogen and 
oxygen in the unperturbed atmosphere. The model is based on 
that of McElroy et al . (1974) , which was modified to give 
self-consistent solutions for a closed system of reactions 
(1) through (7) in photochemical equilibrium above 35 km. 
Reaction rates used are those quoted by McElroy et al . , 
who also supplied photolysis and heating rates, which were 
altered to fit the modified model. 

Fig. 2. Perturbation coefficients for LtO^] and £;,T. 

For i the dashed line is defined by equation (29); the 

so'' id line is *{ given by equation (39). For ‘f the dashed 
line is =* -ilT/Tl, v/here is given by equation (33); 
the solid line is i = i f’ from equation (40). The and 
are independent of 6 except for attenuation of sunlight 
by perturbations overhead. The final ')> and t do depend on 
which was taken uniformly through the atmosphere as 6 = 
+0.1. Below 35 km the t's computed by this technique are 
unreliable, because Dickinson’s (1973) linear perturbation 
treatment of radiative cooling is net approi^riate . Below 
35 km mixing is more rapid than photochemical relaxation and 
we take the t perturbation coefficients as constant and 
fixed by the 35-km values. 
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